Cotunneling renormalization in carbon nanotube quantum dots 



Gediminas Kirsanskas, Jens Paaske, and Karsten Flensberg 
Center for Quantum Devices, Niels Bohr Institute, 
University of Copenhagen, DK-2100 Copenhagen 0, Denmark 
(Dated: March 7, 2013) 
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"shell." It is demonstrated that otherwise degenerate and equally tunnel coupled K and K' states 
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I. INTRODUCTION 

Inelastic cotunneling^ provides a powerful spectro- 
scopic probe of the electronic and magnetic states of 
a Coulomb blockaded nanostructure, be it a quantum 
dotpM' or a single moleculeP^ The basic cotunneling 
process itself corresponds to a virtual charge fluctua- 
tion of the dot, and as such it also ren orma lizes the 
very spectrum which is being investigated!^^] p r suf- 
ficiently strong tunnel couplings this may even instigate 
sharp many-body resonances in the system, such a s the 
celebrated Kondo effect for spinful nano structures] 13 * 14 * 
which in itself is a telling spectroscopic signature. 

For a dot in a specific charge state, the thresholds for 
inelastic cotunneling correspond to energy differences be- 
tween excited states and the ground state. Therefore, 
a threshold is only renormalized by cotunneling if the 
energy of the corresponding excited state is renormal- 
ized differently from the ground state, i.e., if the two 
have different tunnel couplings. Furthermore, since fluc- 
tuations to the two adjacent charge states are generally 
different, this often leads to a marked gate voltage de- 
pendence of t he co tunneling thresholds. This effect was 
first predicte d 15 * 16 * for systems with ferromagnetic leads, 
where charge fluctuations become spin dependent and 
consequently induce a gate dependent exchange field on 
the dot. Experiments have later confirmed the presence 
of this, potentially sizeable, exch ang e fielcP^ as well as 
its characteristic gate dependence}^ and the same effect 
has recently been shown to have a strong influence also 
on the measured tunneling magnetoresistanceP^ 

Experiments with carbon nanotube quantum dots con- 
tacted by either normaP^l or superconducting^ leads 
showed a characteristic gate dependence of inelastic co- 
tunneling thresholds repeated over many different charge 
states. This systematic behavior could again be ascribed 
to cotunneling renormalization, by invoking different tun- 
neling amplitudes for the two different states involved in 
a carbon nanotube "shell." Since, however, K and K' 
states (see below) in a carbon nanotube are time-reversal 
partners, they will have the same Hamiltonian overlap 




Fig. 1. (Color online) Qualitative picture of tunneling in- 
duced JfTf'-mixing for a carbon nanotube coupled to normal 
leads in a zero magnetic B — field, a) Example of the 
tunneling process, which introduces KK'-mixing with a sin- 
gle electron in a carbon nanotube "shell". Here [il,r denotes 
chemical potentials of the leads, b) Schematic of the cotunnel- 
ing threshold (solid and dashed blue lines) for a single charge 
Coulomb diamond, in the absence of spin-orbit coupling and 
disorder induced K K' -mixing, and with real tunneling ampli- 
tudes tK, t K i . The effective tunneling renormalized ground 
state of the quantum dot is \K+) with larger tunneling rate 
Tk+ ~ l*-ff+| 2 (dashed blue line) at the left side of the dia- 
mond and \K—) with smaller tunneling rate Yk- ~ \^k-\ 2 
(solid blue line) at the right side. 



with the lead states and hence equal tunneling ampli- 
tudes. The two differently coupled states observed in 
Ref . therefore must correspond to a mixture of K and 
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K' , reflected also in the non-zero cotunneling threshold. 
This in turn implies that the observed gate dependence 
of thresholds should disappear when applying a strong 
magnetic field along the tube, again making K and K' 
good quantum numbers, which indeed was confirmed in 
a more recent experiment p2l 

In the present paper, we revisit the problem of cotun- 
neling renormalization in a carbon nanotube quantum 
dot, and include now the cotunneling induced mixing of 
the dot states. We also include the effects of spin-orbit 
coupling as well as intrinsic KK'-xnixhxg, A.xk' j due 
to impurities and/or nearby gates or substrate breaking 
the lattice symmetry of the tube. We model the leads 
by a simple flat band metal which carries no informa- 
tion about the quantum numbers, K and K', of the dot 
thus giving rise to tunneling induced KK'-mixing (see 
Fig. [l^) . This assumption excludes from our survey the 
possibility of orbital quasi-degeneracy and the associated 
SU(4) Kondo effect, which has been discussed in connec- 
tion to experiments on certain carbon nanotube quantum 

dotspma 

Using quasi-degenerate perturbation theory, we de- 
termine the leading-order tunneling renormalization of 
thresholds for inelastic cotunneling, strictly valid only 
deep inside the Coulomb diamond delimiting the regions 
of definite dot charge. Throughout the paper, we ne- 
glect the sub-leading effects of level broadening. While 
formally correct, it should be kept in mind that the tun- 
neling induced fine structure which we report here will 
eventually be slightly smeared in any experimental real- 
ization. A very recent worlP^ has reported interesting 
effects of tunneling renormalization of the lead-dot tun- 
neling amplitudes themselves, and also this effect will 
be omitted here as a subleading effect. Figure [IJd il- 
lustrates one of our central findings, namely, the ten- 
dency for otherwise degenerate dot states to split up 
away from the particle-hole symmetric midpoint of the 
Coulomb diamond. Since K and K' states have equal 
tunnel couplings, the system takes advantage of the tun- 
neling induced KK' -mixing allowing to form respectively 
a strongly coupled ground state, \K+), near the bound- 
ary for removing an electron (left side), and a weakly 
coupled ground state, \K— ), when close to adding an 
electron (right side). Altogether, this gives rise to a char- 
acteristic gate voltage dependent splitting of these nom- 
inally degenerate tube states. With outset in this sce- 
nario, we will then include the extra effects of spin-orbit 
coupling, intrinsic if if '-mixing, and an applied magnetic 
field. 

The paper is organized as follows. In Section [IT] a con- 
stant interaction model for a carbon nanotube quantum 
dot is introduced. In Section |III| the tunneling induced 
shifts of eigenenergies are considered perturbatively to 
second order, and in Section [IV] effects on cotunneling 
thresholds due to these energy shifts are discussed. In 
Section [V] we examine the magnetic field dependence of 
the renormalized cotunneling thresholds, and concluding 
remarks are made in SectionlVll A more elaborate discus- 



sion of the single particle states and the tunneling ampli- 
tudes in a carbon nanotube quantum dot is relegated to 
Appendix [X] and comparison of second- and fourth-order 
perturbation theory in the tunneling Hamiltonian is pro- 
vided in Appendix [B] In Appendix [C] detailed derivation 
ofEq. ( 19]) is given. Throughout the paper, except in Ap- 
pendix ~M we employ units in which the reduced Planck 
constant, the elementary charge, and the Bohr magneton 
are equal to one, i.e., H = e = /ig = 1- 



II. MODEL 

The system under consideration is a quantum dot con- 
nected to source and drain electrodes and capacitively 
coupled to a gate electrode. Such a system can be mod- 
eled by the following Hamiltonian 



H = H 



LR 



Hjy + Hi 



where 



LR 



OLVS 



(1) 



(2) 



describes the source and drain electrodes as two reser- 
voirs of noninteracting electrons. Here a — L, R stands 
for the left or right lead, Hl,r — ±V/2 is the chemical 
potential of the leads, which depends on the applied bias 
voltage V, s =f, I denotes the spin of the electron, and 
v are other quantum numbers of electrons. The quan- 
tum dot region is described by and we model it as 
four localized single particle states with the interaction 
between electrons assumed to be constant: 



71=1 



n did n + U(N - N g f 



(3) 



Here N g corresponds to the gate voltage (V g — 2UN g ), 
N = Y^ti=i d\id n is the total number of electrons on the 
dot and 2U denotes the total capacitive charging energy 
of the dot. The single particle spectrum e n for a carbon 
nanotube quantum dot, which we will be interested in, is 
defined in Section Hi A[ 

The coupling between the leads and the dot is de- 
scribed by the following tunneling Hamiltonian 



Hi 



E 



(t n c f d 



(t n Y(ftr ) 



(4) 



where i"„ s is the tunneling amplitude from the dot state 
ri to the lead state avs. These tunneling amplitudes 
are defined in Section |IIB| We will treat the tunneling 
Hamiltonian Ht as a perturbation to Hlr. + Hr>, when 
examining its influence on the spectrum of the dot. 



A. Single particle states and spectrum 

Carbon nanotubes hav e a nearly fourfold degenerate 
energy level structure I 4 I 27 H H1 which is due to intrinsic spin 
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(t, I), and so-called isospin or valley index (if, if'). In 
a lot of experimental cases every four states, "a shell", 
are separated by large energy AE comparable or bigger 
than the charging energy 2llW The fourfold degeneracy 
is broken by spin-orbit coupling and disorder. However, 
if there is no magnetic field, the system still has time- 
reversal symmetry and we are left with a pair of twofold 
degenerate states, called Kramers doublets. We label 
these nearly fourfold degenerate states as 



\Kt),\K' l),\Kl),\K' t), 



(5) 



which we will refer to as the ifif'-basis. Without mag- 
netic field the time-reversal partners (Kramers doublets) 
in the above set are the states 



\Kt)=T\K> |), 
\K\)=T\K> t), 



(6a) 
(6b) 



where T is the time-reversal operator. 

Following Ref. (H we will use the following single parti- 
cle Hamiltonian written in if if '-basis for a single "shell" 
of a carbon nanotube quantum dot: 



|ift) 




|ift) 

/ cos£ 
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lift) 
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\K> t) 
\A KK , \ 





E-i.+i 
\K> t) 

\ 

sin£ 



cos C / 



(7) 



where Akk' is if if '-mixing due to disorder, g s sa 2 is the 
electron's Lande ^-factor, £ is the angle between magnetic 
B field and the tube axis, and energies in the diagonal 
are given by 



A so BcosC 

E T , S = TS—— =F Tfforb 7, , 



(8) 



where r = 1) for K(K'), s = +1(— 1) for spin up 

(down), g or b is the effective orbital Lande ^-factor, and 
the upper (lower) sign corresponds to the conduction (va- 
lence) band. After diagonalizing the Hamiltonian ([7]), we 
find the eigenspectrum e n . In our discussion of tunneling 
induced energy shifts, we will use the energy spectrum for 
the conduction band and assume Aso > 0. 



B. Tunneling amplitudes 

To define the tunneling amplitudes in the if if '-basis, 
we need to specify the lead states. We assume that the 



lead states constitute Kramers doublets, which we will 
denote: 



\avf) = T\aPi), 



(9) 



where a = L, R is a lead index. We also assume that the 
Hamiltonian describing electrons in the leads is real and 
that we can choose arbitrary spin-quantization direction 
(for example, there is no spin-orbit coupling in the leads) . 
In the case of no spin-orbit interaction the states ^ can 
be represented in real space as 



\av t) 



«/( r ) 




{r\av I) 





aw(r) 



(10) 



where a av (r) is a real function of coordinates. Here the 
spin-quantization direction is chosen along the tube axis. 
When there is no magnetic field, we can define the fol- 
lowing tunneling amplitudes 



,K'i _ 



Ki _ 



.K't 



(av^\H tot \K t) =t oa/>1 , 
(av \\H tot \K' l) = t* av ll 
(av l\H tot \K I) = t aUy2 , 
(ais t\H to t\K' t)=t! 



(11) 



u av,2 ' 



where iitot is the total single particle Hamiltonian rep- 
resenting leads, dot region, and tunneling barriers. We 
also used time-reversal symmetry of the states. It can 
be shown (using a real space representation of carbon 
nanotube states described in Refs. .30-35]) that the dif- 
ference between tunneling amplitudes t av ^\ and t aVt 2 ap- 
pears due to spin-orbit coupling and that it can be ne- 
glected, when spin-orbit coupling is small compared to 
curvature induced splitting in the carbon nanotube (see 
Appendix fAl , i.e., 



"OCT/,- 



tr 



tr 



(12) 



However, if the states in the leads cannot be chosen to 
be real (for instance, if there is spin-orbit coupling in the 
leads), then the two Kramers doublets can be coupled 
differently (t aUt i ^ t aVt 2), which we will also consider. 
Note that in general if there is spin-orbit coupling in the 
leads, the lead states have the form 



(r\av t) 



cw(r) 



(r\av I) 



-KM 
<C(r) 



(13) 



where a ct!/ (r), 6 Q „(r) are complex functions of coordinate. 
This additionally may introduce non-vanishing tunneling 
amplitudes 



OLV\, 



\H tot \K\)=t 



t* v \ = (av^\H tot \K> |) = -C,3, 



{pLV J \\H tot \K I) = t avA , 
(av \.\H tot \K' t) = -t 



(14) 
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The inclusion of such terms does not change the qualita- 
tive picture of the results presented in the article, so for 
simplicity we will neglect them. 



In general the tunneling amplitudes (111 also depend 
on parallel magnetic field, however this effect is small, 
when magnetic field is small compared to curvature in- 
duced splitting, and we do not consider it in our calcu- 
lations (see Appendix |A| . We also neglect v dependence 
of the tunneling amplitudes t al/1 / 2 —> t a l / 2 - 



where odd powers of expansion in Ht vanish. We specify 
only the zeroth- and second-order expressions^ 



H 



(0) 



LR , TjD \ r 

mm i n mm) u mm' j 



rr(2) _ 1 ttT ttT 

n mm' — 9 / , n ml rl lm' 



E,, 



El En 



(17a) 



E, 



(17b) 



III. TUNNELING INDUCED LEVEL SHIFTS 

We want to determine how the tunneling renormal- 
izes the energy of the many-body states of Hamil- 
tonian Hlr + H& in a subset \m) € {|2?)|LR)} 
using quasi-degenerate perturbation theory ("Lowdin 



where E are energies of the states. We also have changed 
the subscripts D, LR, T into superscripts for convenience. 

simply gives the energy of the state 



Expression ( 17a 
if m 



17b |, 



\m) it m = m . In order to evaluate expressions 
we assume that tunneling amplitudes do not depenc 
(tea's ~~ ^ an d will perform ^-sums in expression ( 17b ) 



partitioning 



" M33I34I 



Here state 12?) is one of the sixteen using the flat band approximation for the leads spectrum, 



many-body states of the dot Hamiltonian i?rj, an d |LR) 
is the zero temperature ground state of the leads Hamil- 
tonian Hlr- We perform a unitary transformation e lS to 
Hamiltonian ([I]) in such a way that for the transformed 
Hamiltonian 



H 



-iS 



He 



iS 



(15) 



the matrix elements (m\H\l) vanish to the desired order 
in ifp, where |Z) £ {|2?)|LR')} are states with a different 
lead state of 22lr than the ground state (|LR) ^ |LR)). 
The matrix elements of effective Hamiltonian ( 15 ) for \m) 
states take the following form 



H n 



H 



(0) 



H. 



(2) 



(4) 



(16) 



i.e. 



/ t CLUS 



h>as L as 



de.. 



(18) 



D 



where ^ = e ay — fi a G [-D . . . D] corresponds to the flat 
band. Here p as is the density of states at the Fermi en- 
ergy for electrons in the lead a with spin s. We also 
assume that the bandwidth D is much larger than other 
energy scales in the problem. After performing the sum 
over intermediate lead states |LR') and also the i/-sums, 
we get the following expression for the second-order ma- 
trix element between states |2?)|LR) and |2?')|LR): 



H. 



(2) 
W 



as 
nn 



tas(ttr^K\^){^\dl\V') (in 
+ (CrC>l4l^D)(^DK,|2?') (in 



E-d - E^ + fi a 



D 



E-d — E,,p n — [i a 



D 



In 



E v > — E^ + [i a 



In 



D 



D 



Ex>' — E^, D — fj, c 



(19) 



Here the sum runs over all possible many-body states 
IV'd) of the dot Hamiltonian Hn- More detailed deriva- 



tion of Eq. (19) is given in Appendix ICI 



We see that the effective Hamiltonian H w H^ +H {2 ^ 
has a block structure to second order in i?x, i-e., only 
the matrix elements between the same charge dot states 
|2?) = \Nl) and |2?') = \NV) are non-zero, i.e., H m ,N'V = 
Hni,ni'&nn' ■ Here we chose to classify the states by the 
charge number N and quantum number I. After diago- 
nalizing the effective 16 x 16 Hamiltonian H, we find the 
new eigenspectrum E to second order in ifp- In the next 
section, we will examine the cotunneling thresholds in 
the middle of the single charge diamond. We will be in- 



terested in the regime where the charging energy is much 
larger than the single particle spectrum energies e„, bias 
V, and tunneling rates r™ s = 7rp a;s |£™ s | 2 , in order for 
our expansion in Ht to be valid. Note that all energies 
are measured in units of T = J2 a j r «j = n Pl2aj\taj\ 2 , 
where t a j are tunneling amplitudes defined by ( |11[ ) , with 
neglected v dependence, and assuming spin and lead in- 
dependent density of states throughout, i.e., p as = p. In 
Appendix [Cj we also consider fourth-order energy shifts 
for the cases discussed in Sections IV and |Vj in order to 
show the validity of the expansion. For the parameter 
regime discussed in the paper, we confirm that fourth- 
order corrections to the energy shifts are irrelevant near 
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the middle of the diamond. 



IV. TUNNELING RENORMALIZED 
COTUNNELING 

In the Coulomb blockade regime with a particular 
charge state N of the dot, an inelastic cotunneling pro- 
cess occurs when the bias voltage V matches the energy 
difference between two different orbital states of that 
charge state (|Vth| = Em — Emv) and a threshold for 
charge transport from the source to the drain electrode 
through the quantum dot is reached. We consider co- 
tunneling processes involving the ground state Em' and 
some excited state Em- In this section, we are interested 
in how this cotunneling threshold is modified by the tun- 
neling in the middle of the single charge N = 1 diamond 
(V g = 2U + AVg/2). We solve the equation 



\V th \= AE(B, V g , Vlh) = E m - E NV 



(20) 



where E is the tunneling renormalized eigenspectrum. 

First we study the case when B = 0, all single par- 
ticle orbitals are coupled symmetrically to the leads 
Tqj = r/4, and the tunneling rate T is much larger 
than the disorder splitting Akk 1 ~ and spin-orbit cou- 
pling Aso ~ 0, with the resulting tunneling renormalized 
threshold shown in Fig. [2^,. We see that the fourfold de- 
generate "shell" spectrum is split in a gate dependent 
way, and we get the 1 st cotunneling line at zero bias (for 
B 7^ it is visible at finite bias) and one twofold de- 
generate gate dependent line (which splits into the 2 nd 
and 3 rd cotunneling lines for B ^ 0). This splitting ap- 
pears due to effective mixing of K and K' states because 
of tunneling to the leads. We can see this more clearly 
by examining the single N = 1 charge 4x4 block of the 
effective Hamiltonian (19 1, which in this case is 



|l,A-t) \l,K'i) \l,Kl) |l,Jf't) 



H d 








H \ 





H d 


H 







H 


n d 




H D 








n d J 



H = 



(21) 

where denotes the single charge many-body eigen- 
states of Hamiltonian (J3j) , with / being the occupied single 
particle state. The off-diagonal term is 



Ho 



r 

-in 



i + 



AV g +aV 
2U 



1 



AV g +aV 
2U 



r av„ 



7T 2U 



(22) 



and Hd is the diagonal term, which is the same for all 
single charge states. Also here a = +1(— 1) for the left 



(right) lead, and in Eq. (22) we have linearized the log- 
arithms in V 2ij aV ■ The Hamiltonian (21) has the fol- 
lowing two-fold degenerate eigenvalues 



E 



H,± TAV ° 



4,1, 



(23) 



with corresponding eigenstates 
1 
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11,^8) = 



l,Ks)±\l,K's)) 



(24) 



On the left side of the diamond {AV g < 0) the ground 
states are the states \l,K + s) and on the right side 
(AVg> 0) are the states |1, K~sV So qualitatively 
Fig. Vlh represents the result in Fig. [5k. From ( 23 ) eigen- 
spectrum we find that the slope of the 2 nd , 3 rd cotunnel- 
ing lines is given by (dashed (red) line in Fig. (2k) 



S 



(25) 



We see that there is a crossing of cotunneling lines ex- 
actly in the middle of the diamond, i.e., the tunneling 
renormalization effects vanish. If the tunneling couplings 
are different for the left and the right leads, the positive 
and the negative slopes get corrections and instead of the 
above expression ( 25 ) , we obtain for the slope 



S+=± 



r, 



1 ± 



ttU 



AV„ 



(26) 



when the tunneling couplings are real, and where r a = 
r Q i +r Q 2. We see that asymmetry between positive and 
negative bias appears, however, this corresponds only to 
second-order effect in T/U. 

When spin-orbit coupling is included, instead of cotun- 
neling lines crossing in the middle of the diamond, we get 
an anticrossing of size (Fig. [21: 



Ai=2[l 



r 



*SO| 



(27) 



It is not possible to restore a crossing if both Kramers 
doublets defined in Section III Al have the same tunnel- 
ing couplings to the leads (t a ,i ~ t a ^), even if mixing of 
K and K' due to disorder Akk' is included (Fig. [2ji) . 
In this case, for symmetric couplings to the left and 
right leads we can find the eigenspectrum around the 
middle of the diamond by linearizing the logarithms in 
Eq. (19), which yields the following energy difference be- 



tween ground and excited states 



AE =< 



— cos(^)— — AV n 
■nu As 



r 

TrfJ 



A 



so 



cos(0) + isin(0) AV; 




where 



As = 



A 2 



(28) 



(29) 



and the phase <j) is the sum of the ifif' -mixing phase 
(Akk' = I Akk 1 \e lc ^ KK ' ) and the tunneling amplitudes 
phase (t L/RA/2 = 1^*): 



>KK' 



■2& 



(30) 



G 




(c) A so = 0, A KK , = 0.16r (d) A so = o.ir < A KK , = 0.13r 



Fig. 2. (Color online) Tunneling renormalized cotunneling thresholds shown as solid curves (blue) for different K K' -mixing and 
spin-orbit coupling values. The dashed lines (black) depict bare cotunneling thresholds and dashed-dotted lines (gray) show 
bare sequential tunneling thresholds. The long dashed lines (red) in a) show the linearized cotunneling threshold, when there 
is no K K' -mixing and spin-orbit coupling, and the dotted lines (red) in c) and d) depict, respectively, the position of zero bias 
crossing and minimum separation of the cotunneling threshold from zero bias. The values of the other parameters are 5 = 0, 
U — 32r, D — 10 9 r, with all tunneling amplitudes t L / R ^/2 being equal and P^g\t L / R1 / 2 ] = 0. 



From the above expression ( 28 ) we find that an anticross- 
ing appears near the point 



ttU , ( 

1 - — ) A KK , cos(0), 



and its size is 



An = 2 1 - 



ttU 



A| +sin 2 ((/))A 



2 

KK' 



(31) 



(32) 



We note that in this case the middle of the bare dia- 



mond is given by AV„ 



-As- If spin-orbit coupling 



is neglected Aso = 0, but there is KK' -mixing Akk 1 , 
and the total phase <fi is equal to zero, a crossing instead 
of anticrossing appears near the point (31), as shown in 
Fig#. 

By changing the relative phases between the left and 
right couplings, we can reduce the tunneling renormal- 
ization, as shown in Fig. [3j The condition for complete 



reduction of the tunneling renormalization around the 
middle of the diamond is 



thstha + tR,ltR,2 
|2 , i 

,1 

which can be rewritten as 



0, 

\t L ,i\ 2 + \t RA \ 2 = \t L , 2 \ 2 



tR,l 



t R - 



(33a) 
(33b) 

(34) 



where £_L/i?,i/2 are complex numbers, and tp is some ar- 
bitrary phase. 

It is possible to get a crossing if the Kramers dou- 
blets are coupled differently, i.e., t a> i ^ t a< 2, as shown 
in Fig. [4j where we for simplicity consider the case with 
only spin-orbit coupling, but the statement is also true 
when K if '-mixing is included. In this case, the condi- 
tion for crossing is given by relation ( 33a ) , when relation 
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I 



1~ 



I 



l.O 



2U 



Fig. 3. (Color online) Illustration of tunneling renormaliza- 
tion reduction, when there is a phase difference between the 
left and right lead tunneling couplings, and the relations ( |33[ ) 
are satisfied. The values of parameters are Aso = 0.1T, 



A KK , = 0.13r, B — 0, U — 32r, D = io 9 r, 



t/2 



Arg[i 



L,l/2 



= 0. The legend is the same as in Fig 




Fig. 4. (Color online) Illustration of a zero bias crossing, when 
there is only spin-orbit coupling, and the Kramers doublets 
have different tunneling couplings to the leads (only relation 
( 33a I is satisfied) . The values of the parameters are Aso = 



0.16F, A, 



B = 0, U = 32r, D 



1CTT, 



Arg[t iil/2 ] = Arg[t fl i] = 0, kxg[t R2 ] = it. The legend is the 
same as in Fig. [2] 



( 33b ) is not satisfied. Then the position of the crossing 



is given by 



AV„ 



-U 



iY 



(35) 



where Tj = T Lij + T R j, 
amond in this case is AV g — 
observe it experimentally, the value of Eq. ( |35[ ) has to be 
between —2U and 2U. 



and the middle of the bare di- 
Ago|- To be able to 



V. GATE DEPENDENCE OF g-FACTORS 

Now we investigate how the tunneling induced level 
shifts affect the magnetic field dependence of the cotun- 
neling threshold, i.e., we examine gate dependence of g- 
factors. As in the previous section, we restrict our exam- 
ination to the single charge N — 1 diamond. We start 
by examining the case when the spin-orbit coupling and 
KK'-mixing are neglected, the magnetic field is parallel 
to the tube axis, and the couplings to the left and right 
leads are equal. When the gate voltage is exactly in the 
middle of the diamond the bare ^-factors are almost un- 
affected and they are reduced only by a factor (1 — Y/ttU) 



3s/orb ~ (1 — k) ffs/orb, 



(36) 



For a very small T/U ratio, the renormalized threshold 
matches the bare one (dashed (black) curves in Fig. [5^,) . 
For perpendicular magnetic field B± the renormalization 
of g- factors is also small and for 1 st , 2 nd , and 3 rd cotun- 
neling lines, respectively, is given by 



K,g s , (1 - K)g s , g s . 



(37) 



Going away from the middle of the diamond, we find that 
g-factors acquires gate dependence (solid (blue) curves in 
Fig. [5^l), which for the nAV g >> B > case is written 



in Table I, for different transitions. The situation when 
spin-orbit coupling is included is depicted in Fig. [5Jd. We 
see that the tunneling renormalization again acts as gate 
dependent Akk' splitting. The effective ^-factors for 
small magnetic B fields {nAV gi Aso >> B > 0) are 
written in Table I, where the following notation is intro- 
duced: 



1- k A 



so 



K AV a 



-1/2 



g s fa g s (l - k) 

fforb ~ 5orb(l - k) 



1 — K A 



so 



-,-1/2 



k AV„ 



K AVg 

1 - k Aso 



-1/2 



(38a) 
(38b) 
(38c) 



When there is only K K' -mixing, the magnetic field de- 
pendence of the cotunncling threshold does not change 
qualitatively, and the only difference at finite AV g is ef- 
fective enhancement of the KK' -mixing due to tunneling- 
renormalization, as can be seen from Fig. [5]:. The situa- 
tion when both spin-orbit coupling and Xi^T'-mixing are 
included is depicted in Fig. [5ji. 

In an experiment, the cotunncling threshold depen- 
dence on the magnetic field angle C, (with respect to the 
tube axis) also could be measured. The angle dependence 
of the bare cotunncling threshold is shown in Fig. [6^,. 
Again, in the middle of the diamond there is almost no 
renormalization due to tunneling and it matches Fig. [6^,, 



Cotunneling 
line 


K,AVg » B > 0, A so = 


nAV g > A so >> B > 


Bare g-factor 


Renormalized g-factor 


Bare g-factor 


Renormalized ^-factor 


Parallel field E 


II 




9s 


(1 - n)g a 


g S + <7orb 


(1 — K — k)g s + <7orb 


2 nd 


Sorb 


K 9s + 2 K AV g 




(1 - n)g 3 


3 rd 


g S + fforb 


9$ + 2fiAV g 


<?orb 


—kg s + g olh 


Perpendicular fiek 


1 B± 







9s 





Kg s + g s 


2 nd 


9s 


ng B 


9s B 
2A SO 


ng a 


3 rd 


9s 


(1 - K)g a 


9s B 
2A SO 


gs 



Table I. Bare and renormalized y-factors for the carbon nanotube quantum dot, when g OT \, > g s and A KK i = 0. 




0.1 0.1 . ^ t i 0.1 0.1 



(c) Aso = a kk , = o.ier (d) A so = o.ir, a kk , = o.i3r 

Fig. 5. (Color online) Comparison of the dependence on parallel B^ and perpendicular B± magnetic fields of the bare (dashed 
curves (black)) and the tunneling renormalized (solid curves (blue)) cotunneling thresholds for different values of K K '-mixing 
and spin-orbit coupling at a gate voltage V g = 1.2 x 2U away from the middle of the diamond. The values of the other 
parameters are U = 32F, D — 10 9 F, with all tunneling amplitudes t L / R ^/ 2 being equal and Arg[ti/a x/2] — 0. 



and the situation at finite AV g is shown in Fig. [BJd. 
If the quantum dot is coupled to the leads asymmetri- 
cally, then the asymmetry between positive and negative 
bias thresholds acquires angle dependence, which can be 
seen by adding positive and negative bias thresholds in 
Fig. [6b, and the result is shown in Fig. [6fc. 



VI. CONCLUSIONS 



In this paper we have examined tunneling renormaliza- 
tion of the quantum dot cotunneling spectrum by consid- 
ering energy shifts of many-bo dy eige nstates using quasi- 
degenerate perturbation theory^UMI j n tunneling Hamil- 
tonian ifr- The second-order result Eq. ( 19 ) is applicable 



to any quantum dot with an arbitrary number of single 
particle orbitals, when the tunneling rates T are much 
smaller than half of the charging energy (r << U), and 
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V 

r 



(a) Bare cotunneling threshold 




(b) Tunneling renormalized cotunneling threshold 
AV 




(c) Asymmetry between negative and positive bias 
tunneling renormalized cotunneling thresholds 



Fig. 6. (Color online) Angle dependence of the cotunnel- 
ing threshold. The values of the parameters are B = 0.13F 

A so = o.ier, a kk , =o,u = 32r, d 

Arg[t L/flil/2 ] = 0. 



io 9 r, = 7i 

t R,l/2 



gate voltage is far from the charge degeneracy points. 

Using this second-order result we determined the en- 
ergy shifts for the carbon nanotube quantum dot, where a 
fourfold "shell" structure (Section II A I of the single par- 
ticle spectrum was assumed. It was shown that tunneling 
renormalization introduces gate dependent KK'-mhung 
of carbon nanotube orbitals, and this in turn renormal- 
izes g-factors for some cotunneling lines in a gate depen- 
dent way. From the energy shifts the cotunneling spec- 
trum was obtained and we found that for asymmetric 
tunneling couplings to the right and left leads bias asym- 
metry appears, which is a second-order effect in small 
parameter T/U. By measuring the cotunneling thresh- 
old asymmetry (if the coupling to the left and to the right 
lead is different) between positive and negative bias and 
its dependence on magnetic field angle with respect to 
the tube axis, it would be possible to indicate that the 
gate dependence of cotunneling lines appears due to tun- 
neling renormalization and not some other effects, e. g., 
a change of the local electrostatic potential. It was also 
found that the tunneling renormalization can be reduced 
by changing the relative phases between the left and right 
couplings to the leads. 



ACKNOWLEDGMENTS 

We thank Kasper Grove-Rasmussen, Jonas Haupt- 
mann, and Martin Leijnse for valuable discussions. The 
research leading to these results has received fund- 
ing from the European Union Seventh Framework Pro- 
gramme (FP7/2007-2013) under the grant agreement no 
270369 (ELFOS). 



Appendix A: KK'-h&sis states and tunneling 
amplitudes 

In this Appendix, we discuss the explicit form of KK 1 - 
basis states for a single "shell" used in Hamiltonian d7|), 
and corresponding tunneling amplitudes ( 11 1 and (Il4l to 



these states. The effective single particle Hamiltonian for 
the carbon nanotube derived using k • p expansion near 
so-called K and K' points i n the first Brillouin zone is 
given in sublattice a space b> |29 | 3Q | 



H = fw f 



Tsk k c - ikt 
k c + ik t Tsk , 



(Al) 



where 



k c = r(k c + — / 
k t = k t - rfc tjCV , 

VT\ 1 

~J) r' 



( m ~y) 



cv - Tsk so , (A2a) 

(A2b) 

meZ, v = 0,±l, (A2c) 



with t = +1(-1) for K(K'), s = +1(-1) for spin up 
(down), the above value of v depending on the type of 
the nanotube, R denoting its radius, and Vp sa 8.1 x 
10 14 nm/s being Fermi velocity in a nanotube. Also k c 
corresponds to wave number along the circumferential 
direction of the tube, which we will denote c, and kt 
along the tube axis direction, which we will denote t. 
The other parameters are given by 



A 



«so 



c 



k 



ai 2^i?' 
a 2 n f C n cos(3#), 



2hv F R 
cos(36>) 



Ahv F R 2 ' 
sin(30) 



4i? 2 
$>ab 1 = 
$o R 



ai = 0.048 nm, 
a 2 = —0.045 nm, 
j3 = 24 meV x nm 2 , 
C = -0.18 nm, 



7re 



AB 



nR 2 B ll , $, 



(A3a) 

(A3b) 

(A3c) 

(A3d) 
h 

7 

e 

(A3e) 



where the terms k , kgo ar e due to hybridization of so- 
called er — 7r bands induced by spin-orbit coupling, fc CiCV , 
fci jCV are induced by curvature of the tube, and fc$ is due 
to Aharanov-Bohm flux §ab through the cross section 
of the tube. Also 9 g [0, 7r/6] denotes the chiral angle of 
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the tube, Ac ~ 6 meV is the strength of the atomic spin- 
orbit coupling, and a±, 012, ft, C are constants calculated 
using the tight-binding approach.^ If the carbon nan- 
otube has a finite length L and is confined by very sharp 
rectangular potential near the ends, th en th e wavevector 
kt along the tube axis becomes discrete^^ 



-, n G Z. 



The eigenvalues of Hamiltonian ( Al ) are 



E T:Si k c ,k t = hv F ±\f k 2 + k 2 + rsk, 



(A4) 



(A5) 



and eigenstates consistent with k • p expansion in sublat- 
tice a space are 



*r : s,fc c ,fe t (c, t) 



^T,s,k c ,kt 



± 



1 

~ V2^L 

k c — ik t 



e iK T T e i(fc c c+i,() / ZT,s,k a .k t 



>k 2 + k 2 

r = {ccos 9 — t sin 9, csin 9 + t cos 9}. 

(A6) 

where the minus sign in z TjS corresponds to the valence 
band, the plus sign corresponds to the conduction band, 
and K+ (K~) is for the K (K r ) point. These are the 
following states, which correspond to the KK '-basis used 
in Section iHAl 

*+i,+i,fc ,fc 4 -> \K t), (A7a) 
*-i,-i,-fc ,- fc4 -> l^'l), (A7b) 

*+i,-i,fc.,* ( -> |# I), (A7c) 
tt-i.+i,-*.,-*, -> |tf't>, (A7d) 



and when there is no magnetic field they are related as 

* + l, + l,k c> fc t = l,-l,-k c ,-kt> 



(A8) 



>> 



We note that — K = K. In the case when \k c — rk c 
|fcso|j|fc$| we can expand the square root in the (A5| 
eigenenergies to give 

E T , s ,k c ,k t ~ E + st—- =f t , (A9) 



where 
E 



±hv F y(rk c - fc c , cv ) 2 + {rk t - fc tjCV ) 2 , (AlOa) 
/ \ 



Arm = 2hv f 



kn± 



sgn(k c 



rk c ) 



so 



rkt — kt 



Sorb 



= 2 



sgn(fc c 



/ Tkt—k t 
\ i~k c k c 



rk c ) evpR 



/ 



(AlOb) 
(AlOc) 



with fin = eh/2m e being the Bohr magneton. Because 
the term E$ is the same for KK'-h&sis states, after ne- 
glecting it, we arrive at the eigenspectrum (TSl) . 



Now using KK' -basis states (A7) and lead states (13 1, 



which are written in spin s space, we get the following 
tunneling amplitudes 



fTS _ 



giK-r e ir(fc c c+/c t t) 



(a£\H to t\Ts) 

(lie 1 ^) J dvA 6s (r)H tot {v) 

(All) 

rArg [k c - fc c , cv - \{k t ~ fct, cv ) - ts&SO + rk^] , 

(A12) 



A £s = 



A vf , t A uU \ = 1 / a* av b, 

u cxi 




(A13) 



where £ denotes either v f or v \.. In general, the lead 
states also can depend on magnetic field B, however, we 



do not consider such a situation. If \k c 



i(k t 



kt,cv)\ » \skgo — k$\, which is in most cases for car- 
bon nanotubes, we can safely neglect the dependence of 
tunneling amplitudes on spin-orbit coupling (fcso) an d 
parallel magnetic field (fc$) to the tube axis. 



Appendix B: Energy shifts with fourth-order 
corrections 

In this Appendix the necessity to use larger model 
space than that of the single state, which is used for usual 
Rayleigh-Schrodinger perturbation theorypl and appear- 
ance of off-diagonal elements in ( 19 ) is discussed. Also 
to show the region of validity of tl 



ihe perturbative expan- 
sion, we consider the energy shifts with fourth-order cor- 
rections included, for the cases discussed in Sections IV 
and El 

The expression for the fourth-order matrix element of 
the effective Hamiltonian, when projected onto the many- 
body states \m) of Eq . (|3l having a particular lead state 
|LR), is given by (Bl|.™If the projection is performed 
to one particular many-body state |m), then the usual 
expansion of Rayleigh-Schrodinger perturbation theory 
is acquired. In the case, when the single particle level 
spacings A of the quantum dot are much bigger than 
or comparable to the tunneling rates T, the expansion 
becomes invalid, because the first term of (Bl) contains 
terms, which are proportional to TA _1 . This situation 
occurs because it is possible to have intermediate states 
|Z), containing the ground state of the leads |LR). But, if 
it is projected to an extended model space, containing all 
many-body states having |LR) ground state of the leads, 
this situation is resolved. However, such a procedure 
introduces off-diagonal elements, which also appear to 
second order. 

The comparison of the energy differences correspond- 
ing to cotunneling thresholds in Fig. [2] is shown in Fig. [7J 
where solid curves (blue) correspond to energy differences 
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Fig. 7. (Color online) The comparison of energy differences corresponding to cotunneling thresholds in Fig. [5] where solid 
curves (blue) correspond to energy differences with corrections up to second order in Ht, and dotted (green) curves correspond 
to energy differences with corrections up to fourth order in Ht- The values of the other parameters for calculation are B — 0, 
U = 32F, D — 10 9 r, with all tunneling amplitudes t^/jj,i/2 being equal and Axg[ti,/ 11,1/2] = 0. 



with corrections up to second order in iJ T , and dotted 
(green) curves correspond to energy differences with cor- 
rections up to fourth order in Ht. From this figure we see 



that for chosen parameters in our calculation we have a 
wide range of gate voltage V g for which the second-order 
perturbation theory in Ht is valid. 



„(4) _ 1 \ -> T T T tjT 
n mm' — n 2^1 ral n W ' n Vl" ' n l" W 



1,1', I" 



1 



1 



LI 



EttT ttT ttT tjT 
£1 ml- n lm" n m"l' rL Vm' 

4 



\E m — Ei)(E m — Ei')(E m — Ei") (E m i — Ei)(E m > — Ei>)(E m > — E\i, 
8 8 



[E m — Ei)(E m — Ei')(E m " — Ei>) (E rn i — Ei)(E m i — Ei>)(E m » — E{) 



1 



1 



(E m ~ El'){E m " — Ei) \E m — Ei E m n — E v ) (E rn > — E l )(E m n — £)<) \E m > - E v E m n - E, 



1 



[Era" ~ Ei)(E m » — Eii) \E m — Ei E m > — Ey ) (E m — Ei)(E m > — Ei>) \E m » — Ei E m n — Ei 



(Bl) 
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Appendix C: Explicit derivation of Eq. (19 1 



In this Appendix we present more detailed derivatio n of Eg. (19). After inserting the tunneling Hamiltonian Q 
into the second-order effective Hamiltonian expression (17b) and setting m — |2?)|LR), m! — |2?')|LR) we obtain 



H. 



(2) 



\ E E 



|LR')^|LR) nn' 

IV>D> 



E T t-, -f- E,t~. — E T t-. , — E , E T t-> 4- , — E T n , — E , 

LR T> LR/ i/'D LR T>' LR/ Vd 



OL V i 



\ E 



Tin 

W>n> 



(LR|(X>[t2 !/s cL s d„|^D>|LR'><LR / [^D|tS'^ s '^ /i/ , s , ( i n ,p')|LR> 

+ <LR|(©|^ s cL s d„IV'D>|LR'><LR'|(^ D |(^;^ s ,)*<^ QV ' s 'l C ')l LR ) 
+ <LR|(©|(t^J*4^ 3 |^D>|LR'><LR'|^ D |^; iy , s , C t ;i/ , s ,d n ,p'>|LR> 

+(LR|(2J[(tS M )*4c w .|^>[LR0<LR'[(^|(t^ v ,.0*< f c o ,^,,[2J0|IJ^ 

1 1 



*S, s (tl s )*^ s <mi^ D ><fe]4/P'> 



HC„ a yCls(l-n a „s)(V\dUi, u )(i, u \d,\V') 



+ E T>- E i, D 



- ^ - E, 



(CI) 



where n. a!/s = (LR| cj^c^ J LR) € {0, 1} denotes the number of particles in the many-body state |LR) with quantum 
numbers otvs. After the first equality in the above expression the first and the fourth terms in the square brackets 
vanish, and for the second and the third term we have to satisfy a'l/s' = avs, which for the second term gives 
£"lr — -Elr/ = £uu and for the third term gives Elr — ^lr' = — £w We note that |2]) corresponds to the grand- 
canonical Hamiltonian, and the term with the chemical potential fi a is only included when performing thermal averages 
over the lead states. 



Now we will perform the thermal average of the expression (CI) over the lead states using the grand-canonical 
ensemble: 



rr(2) _ w rr(2) 

T>T>' Z^i LR CP',LR 2 



|LR> 



CX 1/ S 

nn 

\4>v) 



tS«(*S«)*/a^K[ite>(itel<,[^> 



'-au + E T>~ E i/j D £ av + E T>I ~ E ^ a 



(C S )**SL(l-/a)<»l4l^D>^D|<i„,P'> 



-z av + E v , - 



(C2) 



where 



fa f\^au Ma) 



»&(e™-n) _|_ I 



(C3) 



is the Fermi-Dirac distribution and 



_ i_„-/3 L (-E L -A'LJVL)„-/3R(-ER-AtRA r R) 



(C4) 



is the probability to be in particular lead state |LR) with /?l r being the inverse temperatures, A/j^R being the number 
of particles, and -E?l,r denoting energies of the left and the right lead, respectively. Also 



2 _ e -/?L(BL-MLW L ) e -/3n(BR- A1 RjV R ) 
|LR) 



(C5) 



denotes the partition function. 

By assuming that the tunneling amplitudes do not depend on the quantum number z/, approximating the lead 
eigenspectrum by a flat band (i.e., £ = e ai/ — fi a G [-D. . . D]) with constant density of states p as , and taking the 
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inverse temperatures of the leads to be equal /3l = Pr, — ft, the expression \G2\ becomes 

f -D 



H, 



(2) 



CKS 

nn 



2 ^— ' 

as 

nn 
l*D> 



*S.(*S.)*/K)<2?KI^XitoML|i''> 



i 



+(*2J**« s {i - /(e)}<»l4l^D><^Di<i n ,p'> 

Be - -E^ D + Ma 



1 1 

*" -£ - Mc + E. 



Cs(C a y(V\d n W D )^ D \dl\T>') (In 



+(0**Sl^l4l^D>{VDl<i„,P'> In 



— i?^, D + /i Q — D 
Ed — E^ D — iL a 



Ex> — E^ D — fi a — D 



+ ln 



-In 



E-jyl — E^j a + [l a 



E-p' ~ E^ D + [i a — D 



E-d' — E^-p — Ha 



E^ u — V-a — D 



(C6) 



where V denotes the principal part of the integral. Here after the second equality we have taken the zero temperature 
limit (j3 —> +oo), which also corresponds to taking |LR) to be the zero temperature ground state of the leads. For 
very large bandwidth compared to other energy scales we acquire Eq. (19). 
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